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Introduction 



The theory of random geometries provides a useful framework to describe 
a wide spectrum of problems in statistical physics ranging from physics of 
polymers, membranes and domain walls, to problems where the dynamical 
geometry arises as a purely mathematical object. Another area of applications 
of the theory is related to the quantization of the geometrical theories like the 
string theory or the general relativity |^, |^, ^, |^, |^, |^ . 

Nowadays, one can not imagine physics without the quantum theory or 
the general relativity. Both theories describe physical phenomena with a re- 
markable accuracy and both have proven to have great predictive power. Each 
of them has its own domain of applicability and so far there is no experiment 
which would contradict either of them. Both theories do not interfere with 
each other. Theoretical consistency requires, however, that there exists a cov- 
ering theory which would unify the quantum theory and general relativity. 
This means that either gravity must be quantized or quantum theory "grav- 
itized" . The formulation of such a theory is one of the greatest challenges of 
theoretical physics. It attracts researchers from different areas of physics and 
results in numerous independent approaches p. 

The theory of random geometries reported here is a generalization of the 
conventional Euclidean path formalism successfully used as a nonperturbative 
method in field theory. The great interest in the theory of random geometries 
in the last decade was triggered by string theory (see for review |10|, 



The representation of the quantum amplitudes for strings in terms of the am- 
plitudes for 2d gravity coupled minimally to matter fields evolved, in parallel 
to the string interpretation, as a theory of 2d quantum gravity [|l], H. The 
success of the two dimensional theory O, [T3|, 111] and especially of the dy- 



namical triangulation approach |T5|, 0, [T^ which, in particular, allowed for 
addressing nonperturbative questions [0, ^ , and for calculating invariant 
correlations ||2l], challenged researchers to generalize the ideas to higher 
dimensional gravity. The generalization was done stepwise : first to three 
dimensions [||, |4|, |5|, H, then to four |, g g^, |2|, ||, |0|. 



In this review paper we focus on the discretized approach combining 
the lattice regularization with the standard concepts of critical phenomena in 
statistical physics. The paper is organized as follows. After a short intro- 
duction where we recall the basic concepts and the discretization scheme, in 
the successive sections we discuss the statistical physics of one dimensional 
simplicial complexes (branched polymers), the theory of random surfaces and 
four dimensional simplicial gravity. 



The model of branched polymers is solvable [31, 32, 33, 34, 35|. It un- 
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dergoes a phase transition [3^| related to the collapse of geometry and to 
the appearance of singular vertices . An analogous phase transition is also 
encountered in models of random surfaces |31|, ^ and higher dimensional sim- 
plicial gravity | 



The mechanism of the transition can be mapped 



onto the condensation of the balls-in-boxes model Q . 

The model of dynamically triangulated surfaces is also well understood. 
It has three phases : collapsed geometries, branched polymers and 2d Liou- 
ville gravity pO| , |41| , ^]. The gravity phase has a well defined continuum limit 
corresponding to the quantum Liouville field theory. The theory is analyti- 
cally solvable by means of the continuum formalism [|1^, |1^ and by the 
discretized approach using matrix model techniques |16, 45, The scaling 



and universal properties of this phase are well established. 

Our understanding of four dimensional simplicial gravity has improved 
recently. We have gained insight into the phase structure of the model and the 
mechanisms governing the behaviour of the system. Nevertheless, we are still 
far from achieving the ultimate goal of the study, namely, the determination 
of the relation of simplicial gravity to the continuum physics. 

The basic difficulty encountered in investigating the model is the lack 
of analytic methods of summing over four dimensional geometries. For the 
time being the only way of studying the model is the Monte Carlo technique 

BUI- 

By means of this method, the basic properties of the model have been 
determined. We discuss the state of art and summarize the main properties 
of the model in the section on four dimensional gravity. We show that the 

J m We discuss 

The model has two phases 



model possesses a well defined thermodynamic limit |^ 
the phase structure 



the collapsed phase with 
infinite Hausdorff dimension and the elongated phase with the Hausdorff di- 
mension equal two. It was believed that the phase transition between the 
elongated and the crumpled phase under a variation of the gravitational cou- 
pling constant was second order. Massive numerical simulations showed that 
the transition is however discontinuous, meaning that one can not associate 
a continuum physics with the critical point pi], |52[. The discontinuity of the 



transition may be explained in terms of the constrained mean field scenario 
^ . A physical explanation advocated in is that the conformal mode 



gets released at the transition due to the entropical dominance of spiky con- 
figurations, similarly as above the c = 1 barrier in two dimensions p]| , |HT| , |52| . 
According to this, if one extends the model by adding matter fields, there 
may exist another phase, like the Liouville phase in two dimensions, free of 
this instability. Indeed, recent numerical investigations of 4d simplicial gravity 



interacting with vector fields support this scenario |^ 
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Apart from the main line of presentation we discuss the balls-in-boxes 
model which serves as a mean field approximation for the dynamical lattice 
models |3^, |39[ . We also sketch ideas underlying the Monte Carlo simulations 
of dynamical lattices. We end the paper with a short summary. 



1 Preliminaries 

One defines the partition function on the ensemble of geometries {G} : 

Z = J2W[G], (1.1) 

G 

where Ty[G'] is a nonnegative weight function given by the Gibbs measure. In 
the quantization procedure of geometrical theories, the weight is 

W[G] = e-^[^] (1.2) 



where S[G] is the Euclidean action. In this case, the statistical sum (|1 . 1| ) 
corresponds to the quantum amplitudes. The simplest example of a model 
belonging to this class is a model of random paths. The sum over random 
paths gives the free particle propagator as one expects from the quantum 
theory [|l|. 

In general, the construction of the theory follows the ideas of statistical 
field theory. There are however many detail differences. A standard field 
theory is defined on a manifold with an inert geometry. There is no such fixed 
underlying structure which would serve as a reference manifold in the models 
of random geometry. Geometrical quantities like a distance or a Hausdorff 
dimension are the dynamical properties of a given ensemble, not of a single 
manifold as in field theory. In field theory one defines field correlators of 
the type (0(x)0(?/)) where x,y are the points on the basis manifold. One 
investigates then the behaviour in terms of the distance \x — y\. From this 
behaviour one can learn about the excitations in the system. The task is 
more difficult in the case of the random geometries where one cannot fix the 
points X, y because manifolds fluctuate. It is even hard to define the correlators 
between the pairs of points at a given distance r, since the distance between 
points is a global property of the manifold. The correlation functions can be 
found analytically only in few particular cases like branched polymers P^, 



or two dimensional pure gravity. In the latter case the calculation requires an 
elaborate technique which was developed only for this purpose [^, . 



There are many new interesting phenomena not present in field theory, 
like the geometrical collapse, the back reaction of matter fields on geometry 
or the change of the dimensionality of the system. 
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The class of geometries which can be considered in the partition function 
( |1 . 1|) is not hmited to paths, surfaces or higher dimensional hyper-surfaces. 
Geometry can be thought of as any set with a given metric structure such as 
a diagram with a function specifying the distance between vertices. 

A theory is said to be geometrical if the action S[G] is a function of 
geometry only. In calculations, one usually uses redundant representations of 
geometry. Redundancy manifests itself as a gauge symmetry of the action. For 
random paths, for example, a path between two points a, b can be represented 
as a continuous map of the unit interval onto a target space where the path 
is embedded : t ^ x[t], a = x[0], b = x[l]. The path will not change if one 
changes the map to : t ^ y[t] = x[f{t)] where / is a monotonic diffeomorphism 
preserving the ends of the unit interval /(0)=0, /(1) = 1. A geometrical 
action will not change either : 5'[a:;[t]] = S'[?/[t]]. This diffeomorphism invariance 
corresponds to the gauge symmetry of the representation. 

In general for any representation, a set of maps which represent the 
same geometry G, defines an equivalence class (gauge orbit). If one represents 
geometry by a metric tensor g^^,, the geometrical action must be invariant with 
respect to diffeomorphism. The simplest invariants yield the Einstein-Hilbert 
action : 

s[g,.] = xj d"^^ - ^ / ^'^VaR (1-3) 

where R is the scalar curvature. The symmetry of this representation divides 
the metric tensors into diffeomorphism classes. The sum (|1 . 1| ) must be defined 
in such a way that the over-counting of metrics from the same diffeomorphism 
class is avoided. In general there are two strategies to do this. Either one 
picks up only one representative of each equivalence class, which is usually 
technically impossible, or one sums over all elements of each equivalence class 
but at the same time one divides out the volume of the equivalence class (gauge 
orbit) : 

Z = Y.W[G] = T.7Lw[T] (1.4) 

G r 1.^ J 

A practical realization of this idea was elaborated in field theory as the gauge 
fixing procedure. The role of the factor 1/C[T] is played by the Fadeev- 
Popov determinant. The gauge fixing procedure was successfully carried out 
for random paths and random surfaces corresponding to 2d gravity interacting 
minimally with conformal matter with c < 1 0, |12|, |13], 

For continuous geometries, the expression ( |1.4| ) has only a symbolic 
meaning. One has to express the sums in ( |1.4| ) in terms of well defined math- 
ematical quantities which would uniquely specify what is really meant by the 
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sums. To this purpose one introduces a short distance cut-off. In the contin- 
uum approach, one develops the theory covariantly and in the end one intro- 
duces the cut-off to compute the final covariant integrals. On the contrary, 
in the lattice regularization, one introduces the cut-off at the very beginning, 
before starting the whole machinery. In this way one breaks symmetries of 
the continuum theory but one hopes to recover them in the continuum limit 
when the cut-off is carefully sent to zero. The lattice, which is an auxiliary 
construct, is eventually removed from the problem in this limit. 

The idea of a discretization is commonly used in the statistical field 
theory. It is also not new in the context of geometrical models where it is 
known under the name Regge calculus [Q. The Regge's idea was to use a 
piecewise linear manifold to approximate the continuous geometry. Regge's 
lattices have a fixed connectivity. Geometrical degrees of freedom are encoded 
in the link lengths. This idea has proven to be very helpful on the classical level 
where it is just the finite element method to solve the classical field equations. 
Its applicability to define the sum ( |1.1D over random geometries is, however, 
limited by the integration measure problem. In the two dimensional case, the 
integration measure can be deduced from the Liouville theory. It turns out 
that it is given by a highly nonlocal expression involving all link variables 
of the triangulation [^. In general, the measure is not known for Regge 



manifolds. Any attempt to mimic the measure by local expressions fails ||66|| . 
In other words the Regge method is well suited as a method to approximate 
Riemannian structures but not to approximate the sum over them. 

An alternative discretized approach to random geometries is based on 



dynamical triangulations [|T5], |T^, |T7|. The randomness of the geometry is 
encoded in the fluctuating connectivity of the lattice. Link lengths are fixed. 
The dynamical triangulations method is not as good as the Regge calculus 
as an approximation of the classical equations but it is very well suited to 
problem of summing over random geometries. 

The basic Ansatz is that the sum over diffeomorphism classes of the con- 
tinuum approach can be regularized as a sum over equilateral triangulations. 
The volume of the symmetry group for such a triangulation with labelled 
vertices is equal to the number of possible relabellings C[T] = N\ : 

If one uses the non labelled triangulations T instead, the symmetry factor C[T] 
is equal to A^! divided by the number of distinct labelling of the triangulation. 
For triangulations without any symmetry, the factor C[T] is equal one. 
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The sum over dynamical triangulations can be done analytically by the 
matrix model technique. The results of the matrix model and the Liouville the- 
ory agree. This is usually treated as a strong indication that the two methods 
provide correct definitions of the integration measure over random surfaces. 
An advantage of the lattice method is that without changing the Ansatz ( |1.5| ) 
one can generalize it beyond the Liouville phase of two dimensional gravity, in 
particular, to higher dimensional gravity. 

If one applies the Regge calculus to the discretization of the Einstein- 
Hilbert action ( |1.3| ) on the equilateral simplicial lattice one obtains : 

S = KdNd-Kd.2Nd-2 (1.6) 

where and Nfi-2 denote the number of cZ-simplices and {d — 2)-simplices of 
the simplicial manifold. The coupling constants k's are related to the couplings 
of the continuum action ( |1.3|) and to the lattice spacing in the naive continuum 
limit. 

One can extend this discretization procedure to geometrical actions with 
higher derivative terms p7| , or actions which describe interaction of gravity 



with some matter fields |T0|, |7l| . This discretization scheme leaves some 

freedom because one can add to the discretized action terms, which disappear 
in the naive continuum limit. One should perhaps look at the discrete models 
from a different perspective, in which the model is not viewed as a discretiza- 
tion of a continuum theory but rather as a primary definition of the theory. 
Then if one finds a continuum limit then one can ask what is the underlying 
continuum theory related to this limit, and whether this theory is related to 
gravity. In fact this is the most difficult part of the procedure, since we do 
not have the experimental data. All we know about the continuum theory is 
that it should reproduce general relativity in the classical limit. This, together 
with the self consistency requirement, provides the only available checks for 
the consistency of the constructed theory. 



2 Branched polymers 



The model of branched polymers provides a useful ground for testing various 
ideas. It is simply and solvable. It plays a similar role for random geometries 
as the Ising model for statistical field theory ||31], |3^, The model 

undergoes a phase transition related to the change of the Hausdorff dimension 
and to the collapse of geometry p4|, |35|, |36| . A similar transition is also present 



in the models of random surfaces and four dimensional gravity fS^, ^ 
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Branched polymer is a graph (one dimensional simplicial complex) with- 
out loops. Geometry on a branched polymer is given by the geodesic distance 
between vertices of the graph. The distance is defined as the number of links 
of the shortest path between the points. On a tree-like graph there is only one 
path joining any two points. 

In the simplest case, one considers polymers whose vertices are indepen- 
dent in the sense that there is no direct correlation between the branching 
orders. A branching order is defined as the number of links emerging from the 
vertex. The action for such a branched polymers is 

S[T]=f^N[T]-Y.s{q^) (2.1) 

where N[T] is the number of vertices on a branched polymer T, /i is the 
chemical potential and s{q) is a one- vertex action depending only on the vertex 
order q. The sum in the second term of (|2.1| ) runs over all vertices of the 
branched polymer. The grand canonical partition function defined on the 
ensemble of trees of arbitrary size, reads : 

^(-") = E7^exp(^.(g,)-/iiV[T]) =^^(iV)e-'^^. (2.2) 

T i-' J ieT N 

and can be treated as a discrete Laplace transform of the canonical partition 
function z{N), being the statistical sum over the ensemble of trees with N 
vertices. 

It can be shown that the coefficients z{N) grow exponentially 

ziN) ~ iV^"3 exp(/i,,iV) , (2.3) 

for large N. This means that the thermodynamic limit, ^ oo, is well 
defined. The quantity ficr is the critical value of the chemical potential at which 
the grand-canonical function has a singularity. The power-like corrections 
A^"^"^ determine the type of the singularity of the grand canonical partition 
function for Afi = fi — ficr 0"^ : Z{fi) ~ A/i^""^. The singular part of the 
susceptibility defined as the second derivative of Z, behaves as 

Xifi) = Z"(/i) ~ A/i-^ (2.4) 

at the critical fi. The exponent 7 is called a susceptibility exponent or entropy 
exponent. The exponent 7 is generally used to determine the universality class 



of models of random geometries ||12|, [131, 14 
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The susceptibility is related to the puncture-puncture correlation func- 
tion defined as |E2l 021 : 



Gir. /^) = E 7^^"''^^' E ^(^ - dia, b)) (2.5) 

T ^i-'-i a,b 

The internal sum runs over all pairs of points. The delta function selects 
contributions from pairs at a distance r. Integrated over r, the two-point 
correlator gives : 

x(/i) = EG(^/^) (2-6) 

r 

For large r, the correlation function G{r, fi) falls off exponentially and this is 
a general feature of models of random geometries ||7^. One associates a mass 
with the exponential fall-off : 

m = - lim In . (2.7) 

r— >oo f 

This mass is directly related to the large limit. It is physically important 
if the mass scales to zero when Afj, = // — /icr — >■ 0"*". If it does, the geometry 
has a well defined Hausdorff dimension. The dimension is related to the mass 



critical index given by the scaling formula |T2l| : 

m~A//'^«. (2.8) 

If it does not scale, the geometry is collapsed, as we show later. To see that 
the exponent dn may be indeed identified with the Hausdorff dimension in the 
scaling case, it is convenient to consider a counterpart of the puncture-puncture 



correlator ( |2.5|) in the canonical ensemble with a fixed size N. Similarly to 
the partition functions ( |2.2| ) one can relate the correlation functions by the 
discrete Laplace transform : 

G(r,/i) = 5:G(r,iV)e-'^^ (2.9) 

N 

where G{r, N) is the correlation function in the canonical ensemble. Defined 
this way, G{r,N) has an unnatural normalization proportional to z{N). One 
can get rid of it, by defining the normalized correlator : 

= = ^ ~ ^^^^^ ^^-^^^ 
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where the averaging {. . .)n is over the ensemble of trees of size N. Now we 
are interested in the large N behaviour of the function g. This behaviour can 
be extracted from the inverse Laplace transform of G{r,n) ~ exp(— A/^^/'^-f^r), 
( p.5|) for A/i — 0. Small A/x corresponds to large ~ l/^A'- the transform, 
whence the normalized correlation function g{r, N) must be a function of the 
argument r/N^^'^" for large N. The function g{r,N) measures the average 
number of points at a distance r from a random point. Summed over r, 
g{r, N) gives the number of all points N. If one inserts the universal argument 
r/N^/'^" one obtains, that the average distance between points : 



behaves for large as 



{r) = ^J2^9ir,N) (2.11) 



(r) ~ N^/'^" . (2.12) 



The last formula relates the size of the system to the typical linear extension 
of the system, naturally leading to the interpretation of the mass exponent dn 
( p.8|) as the Hausdorff dimension. 

Since the model is solvable, the critical indices : the susceptibility expo- 
nent 7 and the Hausdorff dimension dn can be calculated. 

In practical calculations one considers the ensemble of planar rooted 
trees |3^. Rooted trees have one marked vertex with one attached link. The 



planarity means that trees are drawn on a plane. The existence of the root 
and the planarity uniquely specify the symmetry factor C[T]. One can find a 
representation where C[T] = 1. 

The grand-canonical partition function (^.2|) is given by the solution of 



the equation |31, 34 



^, = K{Z) = log J2piq)Z'^-' (2.13) 

q=l 

which is to be solved for Z. The series coefficients p{q) are related to the 
one- vertex action : p(q) = e~^^'^\ A vertex with order q contributes p{q) to the 
total weight of the tree. The weights p's are nonnegative in unitary models. 
Vertices with the order q are forbidden if p{q) = 0. One requires p{l) > 
in order to have the endpoints in the tree. In order that the polymer could 
branch, at least one weight must be positive for g > 3. 

The solution of the equation (|2.13|) for Z is given by the inverse function 



of : Z = K ^(yu). In fact, without inverting one can find the singularities 
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of Z. Namely Z is singular at a certain /i^ when either the derivative of the 
inverse function is zero, K' = 0, or the inverse function K is itself singular at 

Z,r = Z{^I,r) 11. 

In the former case one gets : 

= K{Z,r) + ]^K"{Z,r){Z - Z,,f + ... (2.14) 

which gives Z ~ A/i^/^ where Afi = fi - K{Z„). It follows that 7 = 1/2, 
as seen from the comparison with the singularity Z ~ of the partition 

function^. This is a generic situation since the function K has a minimum for 
any choice of weights p's fulfilling the general assumptions discussed before. 
The function K is singular when the series ( |2.1tj| ) has a finite radius 



of convergence. For example, for the weights, which for large q behave as 
p{q) ~ , the series has a singularity (1 — ZY~^ at Z = 1. When P > 2, the 
situation is interesting since there are three different possible behaviours of the 
singular part of the grand canonical partition function. Let Zq be the value 
at which the derivative of K vanishes : K'{Zo) = 0. Now we have two special 
points Z = Zq and Z = 1. When Zq < 1, the singularity of Z comes from 
inverting K around Zq. In this case one obtains the generic value 7 = 1/2, 
as previously. When Zq > 1 the singularity of Z comes from the singularity 
of K at Z = 1. In this case Z ~ A/i^~^ and hence 7 = 2 — /?. This is a 
semi-generic situation. Finally, in the marginal situation when Zq = 1, the 
singularity of the partition function is Z ~ A/ii/(^~i) when 2 < p < 3 or 
Z ~ AyU^/^ otherwise. Hence the exponent 7 is equal to (/? — 2)/(/9 — 1) or 1/2 
respectively. 

It is worth noting that the situation will not change if we add an expo- 
nential prefactor to the weights p{q) : q^'^ e^^'^q~^ . Namely, due to the 
Euler relation, the sum of the vertex (branching) orders over the non-rooted 
vertices is : 

^g, = 2iV-i, (2.15) 

i 

and therefore the exponential term can be absorbed in the cosmological term 
e~'^^ ( p.2|) by a redefinition of the chemical potential /i — /i + 2n. After the 
redefinition the weights have again the original form p{q) ~ q~^ and we obtain 
the previous situation. 



^In a rooted ensemble, like in the one considered here, one vertex is fixed by the root. 
Therefore to obtain susceptibility x needs to differentiate Z with respect to ^ only once 



and not twice as in (2.4) where Z is partition function for a non-rooted ensemble. 
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The generic phase of branched polymers has many universal properties. 
In particular the normalized puncture-puncture correlation function has for 
large N the form : 



g{r,N) = VNagi^) (2.16) 



where g is the universal function given by the formula : 

g{x) = 2x6-''' (2.17) 

of the universal argument x = ar/N^^"^. The form of the function g does not 
depend on the choice of p's as long as the system is in the generic phase. The 
universality tells us that the local properties like the branching distribution 
do not affect the long range behaviour. The change of the weights p's can be 
compensated by the change of only one parameter a in the correlation function. 

The Hausdorff dimension is d/^ = 2, as follows from the universal scaling 
X = ar/N^^"^. Finite size calculations show that the scaling is weakly broken 
for finite N and one should use a shifted argument x = a(r + 5)/N^/'^ ^ 



In the marginal situation the form of the correlation function changes 
and so does the universal parameter [^] . In this case the Hausdorff dimension 
is du = 1/7 which is dn = {(3 - l)/(/3 - 2) for 2 < /5 < 3 or dn = 2 
otherwise. In the semi-generic situation the normalized correlation function 
acquires a mass term e"™^ with a non-vanishing mass in front of the scaling 
piece. The mass m does not depend on N. This means that the average 
distance (|2.12|) is of the order 1/m : 

(r) ~ — ~ const (2.18) 
m 

and does not scale with A^. This may be interpreted as an infinite Hausdorff 
dimension. This phase is called a collapsed phase since it is dominated by 
short branched polymers which do not grow, contrary to the generic situation 
dominated by the elongated polymers. 

The collapse of the geometry is a result of the appearance of singular 
vertices on branched polymers |3^, |3^, |3^. A singular vertex is a vertex with 



an order which grows extensively with A^. The mechanism of the appearance 
of the singular vertices can be described in terms of the balls-in-boxes model 
discussed in the next section. The results of the discussion are summarized in 
the table [l|. 

To end this section let us briefly discuss the topological aspect of the 



model One may extend the class of graphs to the ensemble of graphs 
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phase 


7 




generic 


1/2 


2 


marginal (/3 < 3) 

(2 < /? < 3) 


1/2 

(/?-2)/(/?-l) 


2 

(/3-l)/(/3-2) 


collapsed(/3 > 2) 


2-/? 


oo 



Table 1: Critical exponents 7 and dn for the three phases of the branched 
polymers model. 



with loops. Such graphs can be generated by the perturbative expansion of 
a zero dimensional field theory with a potential containing terms 0*^. The 
coefficients in front of the terms are related to the weights p{q). The tree 
diagrams discussed so far come from the leading term in the loop expansion 
corresponding to the classical tree level. The number of diagrams with an 
arbitrary number of loops is not exponentially bounded so that the entropy 
is not an extensive quantity and there is no thermodynamic limit. There are 
some ideas how to cure the problem as discussed below |TB|, ^ 

If one wants to sum over topologies in (|2.2|) one has to consider topo- 



logical terms in the action. The number of loops L is the simplest topological 
term for graphs. The partition function can be written as : 



L 



(2.19) 



where Zl is the partition for the sub-ensemble with L loops, and ^ is a coupling 
constant for the topological term. The nonexistence of the exponential bound 
for the number of all diagrams implies that the series is not summable. It is not 
even Borel summable. The radius of convergence is zero and therefore there 
are functions which can be added to Z without changing the coefficients Zl of 
the series. Such functions are called nonperturbative modes. In principle the 
number of nonperturbative modes is infinite. The idea is to reduce it as much 
as possible. This is done by the loop expansion of the scalar field theory gener- 
ating the asymptotic series ( |2.19| ). This series incorporates contributions from 
diagrams with any number of loops. One can show that in the double scaling 
limit ^ and t = Afx^^"^ /h = const, the susceptibility is given by the 
Riccati equation in the scaling argument t. This equation uniquely determines 
all coefficients Z^ of the series ( |2.19| ) and has only one non-perturbative pa- 
rameter. Thus the number of perturbative modes gets reduced to only one and 
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the goal of summing over topologies gets partially achieved. Unfortunately, so 
far no physical principle to fix this remaining free parameter is known. The 
existence of the double scaling follows from the fact that the susceptibility 
exponent 7 grows linearly with the Euler number. This was first discovered 
in two dimensional gravity |T8|, [19|, |20[ . One also finds a linear dependence in 
the marginal and semi-generic phases of branched polymers |3H . 



3 Balls-in-boxes model 

In this section we discuss a model of weighted integer partitions - mean field 
approximation for models of dynamical lattices. The model undergoes a phase 
transition which has many common features with the Bose-Einstein conden- 
sation. The integer partitions of the model correspond to the partitions of 
vertex orders of djTiamical lattices. The phase transition relies on a conden- 
sation which favours partitions with one integer proportional to the sum of all 



integers in the partition |36|. This integer corresponds to the singular vertex 



order on the random lattice 37, 38 



The partition function of the balls-in-boxes model : 

Z{M,N)= Pi(li)---Pi(lM)S{qi + --- + qM-N) {3.1] 



describes weighted partitions of balls in M boxes. The function 5() is the 
Kronecker delta. The weight is a product of one-box weights p{qi) which means 
that the numbers of balls q in any two different boxes are independent of each 
other. The independence is weakly broken by the constraint on the total sum 
which prevents the factorization. This constraint makes the model nontrivial. 
For the convenience we assume that each box contains at least one ball q>l. 
If we additionally choose N = 2M — 1 then the partition function (|3.1|) is equal 



to the partition function (|2.2| ) of the branched polymer |]36| , M - The numbers 
qi of ( |3.1| ) correspond to the orders of vertices of the branched polymer. The 
Euler relation X)i qi = 2M — 1 introduces the constraint. 

In the large M limit the partition function of the model takes the form : 

Z{M, N) = e*^^(^)+- (3.2) 

where p = M/N is the average density of balls per box, and /(p) is the free 
energy density per box. The model has a phase transition at a certain critical 
density pcr where the free energy is singular. The value of pcr depends on the 
choice of weights p's. In particular p^r may be moved away to infinity or to one. 
In either case the model would only have one phase. Here we are interested in 
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the situation where pcr is finite and the model has two phases depending on 
whether p is lager or smaller than p^r- For example, for the weights : 



the model undergoes the phase transition at 

„ C(/5-i) 



(3.3) 



(3.4) 



where C, is the Riemann Zeta function. To fix attention and not to make 



formulas too abstract we will keep in this section this particular form (|3.3| ) of 
weights, but the discussion can be naturally generalized to other forms |^ p5| . 
The critical density per is finite for j3 larger than two. When (3 goes to infinity 
the critical density approaches one. When (3 goes to two the critical density 
goes to infinity and eventually the transition disappears when {3 becomes equal 
or smaller than two. 

An alternative way of triggering the transition is to change (3 for fixed p 
as for instance for branched polymers where density is fixed p = 2. 

The free energy density has the following singularity at the phase tran- 
sition : 



dpf ~ < 



for p>3 



(3.5) 



There are logarithmic corrections in Ap for integer f5. The free energy has the 
same type of singularity in the parameter A/5. 

Analogously to branched polymers (see previous section), an additional 
exponential factor in the weights (|3.3| ) p{q) = e'^'^q~^ does not affect the phase 
structure ||36|| . 

It is convenient to consider the dressed one-box probability 7r(g) to see 
what happens in the system at the transition. The dressed probability 7r(g) is 
a probability that a particular box has q balls. One obtains for large M : 



7r(g) 



q-Pf.-M 



for p < Pc 



(3.6) 



+ jT5{q - M{p - Per)) for p>pc 
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where the normahzation factor Af{fi) is 



A^(/^) = E ^"^e-'^' (3-7) 

9=1 

and /i is a positive function of p which vanishes at the transition. One can 
check that the average number of balls per box is indeed p : 

oo 

E = P (3-8) 

q=l 



for the one-box probability 7r(g) given by ( p.6|) . The interpretation of the result 
( p.6| ) is following. In the low density phase (p < pcr), the typical fluctuation 
of the box occupation number q is of the order At the transition, p, 

vanishes, so the fluctuations must be of the order of the system size. Indeed, 
at the transition one box captures a number of balls which grows with M 
as shows the argument of the delta function in the second term of the high 
density formula ( |3.6| ), and the occupation of this singular box gives rise to 
large fluctuations. At the transition the probability vr has the critical form : 

vre.(g) = |^. (3.9) 

Above the transition this form is frozen but it is supplemented by the anoma- 
lous term : 

anomaly = ^5 (^g - M(p - Per)) • (3.10) 

This term is anomalous in the sense that it disappears from the probability 
distribution irlq) if one takes the point-wise limit M — > cx) for each fixed q. 
In this limit, only iTcriq) survives. If one calculated the average ( |3.8|) for such 
a limiting probability distribution TTcriq), one would obtain the wrong value 
Per instead of the correct one p. This means that the anomalous term can 



not be neglected in calculating the average ( |3.8D . The anomaly introduces an 
additional probability 1/M of picking one out of M boxes with M{p—pcr) balls. 
This is the singular box (singular vertex). As we see the anomaly corresponds 
to the condensation of balls in one box which just takes over the surplus of 
balls holding the rest of the system critical. The condensation is similar to the 
Bose-Einstein condensation. The difference between the two condensations is 
that in the Bose-Einstein condensation, particles go to the lowest energy state, 
while here the box must be chosen by the symmetry breaking since the boxes 
are indistinguishable. 
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Figure 1: Evolution of the shape of the dressed one-box probabihty 7r(g) with 
density p. The three curves correspond to densities below, at and above the 
transition. (The model with the weights for /3 = 4) 



The transition to the condensed phase is visualized for the finite size 
system M in the figure |l|. The finite size calculations have been done by an 



improved version of the recursive technique described in ||3g . 

There are some secondary finite size effects to the formula ( p.6| ) as for 
instance that the peak at M{p — p^r) is smeared for finite M or that one has 
to go to sufficiently large M to see the peak depart from the remaining part 
of the distribution. However, the basic features of the solution ( p.6|) that the 
position of the peak moves linearly with M and that its height decreases as 
1/M are already seen for moderate sizes M (see figure 0). 



It is also interesting to consider ensembles with varying density p |^3|, ^ |. 



There are two natural candidates. The ensemble with a variable number of 
balls with the partition function : 

Z(M, p) = Y. Z{M, N)e-''^ (3.11) 

N 
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Figure 2: The dressed one-box probability 7r(g) in the condensed phase for 
p = 2, /? = 4, iV = 128, 512, 2048. 



or the ensemble with the variable number of boxes : 

Z{k, N)=Y, ZiM, N)e^^ (3.12) 

M 

If the sum over N in ( p.ll|) extends to infinity, the problem factorizes to M 



copies of the urn-model [^, |7J]. It does not, however, if there is an upper 
limit Nmax for N. In this case the phase structure is basically the same as in 
the Z{M, N) ensemble. For large n the system is in the low density phase. At 
some critical value of fi, the system enters the phase where one of the boxes 
captures the surplus of balls to maximize N. Analogously in the ensemble 
( p.l2| ), for K above a critical value, the system is in the high density phase 
realized by minimizing the number of boxes to the smallest available number 

If one considers the (k, N) ensemble, it is more convenient to use the 
quantity : 

(M) _ 1 dZjK^N) 
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Figure 3: The distribution of r = M/N in the (k, A^) ensemble in the pseudo- 
critical region for two different volumes : (-0.32184,512), (-0.31910,1024) 
for the model with the weights p{q) = q^^ for (3 = 2.5. 



instead of the balls density. The values of r are limited to the range : Mmm/N < 
r < Mmax/N. The lower limit may be naturally chosen : M^m = 1, to have 
at least one box, and the upper one : Mmax = N, which corresponds to one 
ball per box. When k goes from large negative to large positive values, r goes 
from the lower to the upper limit. In the large N limit the function r = r{K) 
can be found. For k < Kcr the system stays at the lower limit r = 0. At 
the critical value it jumps to ro > and then approaches continuously the 
upper limit when k goes to infinity. The phase transition is discontinuous and 
there is a latent heat related to the height ro of the jump. For finite N the 
transition between the phases is smoothed. There is a crossover interval of the 
size 5k ~ A^~^ where the curve r steeply goes between the two regimes. In this 
crossover region system is effectively a mixture of two phases and hence the 
distribution of r has two peaks as shown in figure ^ One peak is at the lower 
kinematic limit while the other at rg in the large r phase. When k moves in the 
crossover region, the relative peak heights change. One can define a pseudo 
critical value of k for a finite system as the value at which the heights of both 
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peaks become equal. We show in figure ^ histograms obtained by finite size 
computations for two A^'s. The depth of the valley between the peaks increases 
with the size since each of the peaks becomes narrower. Eventually in the limit 
N ^ oo the configurations from the valley are completely suppressed. 

To summarize, the transition related to the appearance of the surplus 
anomaly is continuous in the fixed density ensemble and discontinuous in the 
(M, fi) or (k, A^) ensembles with fiuctuating density. 



4 Random surfaces 

A theory of random surfaces has been an active research field since Polyakov 
proposed the geometrical approach to the quantization of strings by combining 
the Feynman quantization principle with the geometrical nature of the string 
action ^ 0]. There are some excellent reviews summarizing the models, the 
ideas and the methods P, [TO, |TT|. Here, for completeness, we briefiy sketch 



the ideas and the main results which can be found there. Then we discuss 
some issues which appeared later in the literature. 

Let us come to the origin. The partition function for the Polyakov theory 
reads : 

^ = / ^f^^'^ [ ~ ^^9f.u, V?]] (4.1) 

The action is a sum of the Einstein-Hilbert action ( p..3|) and the action for the 
matter fields (p coupled minimally to gravity. The term : / d^C,y/gR = 47r(l — /i) 
of the Einstein-Hilbert action is a topological invariant and appears only if 
one sums over genera h. 

The Nambu-Goto string embedded in D dimensions rewritten in terms 
of the Polyakov formalism corresponds to the action : 

S[gab, ^] = E / d'^^g'^'da^'^d,^'^ + l^fd'C (4.2) 

fi=i 

for D scalar fields coupled minimally to gravity. The term proportional to the 
Euler characteristic is not displayed. One extends this action to non integer 
D. In this case D = c where c is the central charge of the conformal matter 
coupled minimally to gravity. 

As mentioned, the continuum Liouville field theory and the discretized 
dynamical triangulation approach used to calculate the partition function, 
yield the same results for c < 1. The approaches are independent. Even more, 
they are independent to such a degree that it is very difficult to find a direct 
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correspondence between them. It is much easier to compare results for the 
universal quantities than to compare the formalisms themselves. For example, 
a lot of work has been spent to recover the complex structure or moduli spaces 
in the dynamical triangulation approach |^ . 



The fundamental results of the theory are summarized in the KPZ (Kni- 



zhnik, Polyakov, Zamolodchikov) formula [|T^, 13, 14| for the scaling dimensions 
of operators dressed by gravity and the susceptibility exponent 7 defined by 
the analogous formula as ( |2.3| ). Conformal weight Aq of an operator acquires 
a new value A when the operator is coupled to gravity. The dressed value is 
given by the equation : 

A-Ao = -yA(A-l) (4.3) 

where 



« = ^(V25^-v/r^). (4.4) 
The number of surfaces with area A is given by 

where the entropy exponent 7 for spherical surfaces is : 

7 = ^(c-l-v/(25-c)(l-c)). (4.6) 

For other topologies the exponent 7 changes linearly with genus : 7h = 7 + 
h{2 — 7). The Liouville theory breaks down at c = 1. This is known as the 
c = 1 barrier. This is related to the instability of the conformal mode which 
drives the system to the branched polymer phase. In the language of strings 
it corresponds to a tachyonic state which destabilizes the stringy vacuum. 

In the discretized approach, the functional integral over surfaces is reg- 
ularized by the sum over triangulations |T^, [16|, [l^. The discretized theory is 
given by the partition function : 

Z = T.I[<ir7^e-^^^''^exp [-If: E(^r - XfY] (4.7) 



C[T] 



where the sum runs over triangulation, and are D scalar fields with the 
nearest neighbours interactions. For non integer D one can either directly 
weight triangulations by the power —D/2 of the determinant of the Laplacian 
obtained by integration of the X fields or one can consider various statistical 
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Figure 4: The phase structure of the model given by the partition function 
m in the (a, D) plane. 

models corresponding to the conformal matter field with D = c. The number 
of triangles on the lattice (area) is denoted by A, and the genus by h. The 
two terms —fiA + Xh correspond to the Einstein-Hilbert action. 

The product of the powers of vertex orders was originally introduced to 
investigate the stability of the discretized integration measure ^ Q . The 
measure term corresponds to the higher derivative terms and is irrelevant in 
the perturbative regime. 

The phase structure of the model in the {a,D = c) plane has three 
phases : the gravitational phase corresponding to the Liouville theory, the 
collapsed phase with singular geometries and the branched polymer phase 



l|, ^ . The phase structure is approximately sketched in figure ^. 

For c between zero and one, in the Liouville phase, there exists a dis- 
crete series of models of unitary conformal matter coupled to gravity with 
the conformal charge c = 1 — Q/m{m + 1) which can be enumerated by an 
index m = 2,3, .. . @, |Tl|. The susceptibility exponent for this series is 



7 = —1/m ( |4.6|) . This series has a realization in terms of statistical models on 



a dynamical triangulation. For example, the first member of the series corre- 
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sponds to the dynamical triangulation without dressing (pure gravity). The 
second to the critical Ising spins on a dynamical triangulation, the third to 
the three-state Potts model. To associate with a statistical model a particu- 
lar conformal field, one has to compare the operator contents. The conformal 
weights of the underlying conformal theory are related to critical indices of the 
corresponding statistical models. For example for the Ising model on random 
lattice the values of the standard critical exponents |^ |7B|, |77| : 



a = -l,/3= 1/2,7 = 2,5 = 5,^^1^ = 3. (4.8) 



correspond to the exponents calculated from the conformal weights (|4.3|) of 
the c = 1/2 conformal field dressed by the Liouville field. The exponent v 
appears in the combination with the fractal dimension du-, which itself is a 
dynamical quantity. It will be discussed below. 

Analytic calculations for the discretized random lattice are performed 
by the matrix model technique (see reviews 0, |lll)- Let us briefiy recall 
the idea. By the duality transformation one can rewrite the sum over trian- 
gulations as a sum over 0^ Feynman diagrams generated by the perturbation 
expansion of the 0^ matrix field theory in zero dimensions. Amongst diagrams 
of the 0^ theory there are such which include tadpoles and self-energy sub- 
diagrams. These correspond to pathological triangulations containing for ex- 
ample triangles whose two edges are glued together. One can remove tadpoles 
and self-energy sub-diagrams by using standard renormalization procedure for 
the perturbation theory. The renormalized theory has the same universal con- 
tent encoded in the critical indices as the original theory. In the Ising model. 



one can check this by direct calculations DT^j. Universal properties do not 
change when instead of 0^ one considers the 0^ diagrams, as expected on the 
grounds of the more general argument that the local properties of the lattice 
do not matter in the limit when the lattice spacing goes to zero. This intu- 
itively means that in this limit one cannot distinguish whether the lattice is 
built from triangles, quadrangles or other polygons. 

The perturbation expansion of the matrix models generates all terms in 
the partition function ( [4.7|) . The symmetry factor \jC occurs automatically 
from the Wick theorem. The area term corresponds to the perturbation order 
of the diagram which counts the number of vertices. The topological term 



arises from the colour expansion of the matrix field |T^. The matter content 
of the theory is generated by the multi matrix action with chain interactions 



45| , In this way one can construct c < 1 conformal matter field from the 



^Compare with the Onsager exponents for the fixed lattice : a = 0, /3 = f/8, 7 = 7/4, 
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unitary series c = 1 — 6/m(m + l) or some non-unitary matter. Within the for- 
mahsm one can calculate critical exponents, correlation functions, macroscopic 
loop amplitudes etc. The colour expansion of the matrix model simultaneously 
incorporates contributions from all topologies. This leads to a partial solution 
of the problem of summing over genera by reducing the number of nonpertur- 
bative modes similarly as discussed in the section on branched polymers. For 
c = 0, for example, the sum over topologies is reduced in the double scaling 
limit to solutions of the Painleve II equation which has only two nonpertur- 
bative modes [0, [T^, ^ . 



Contrary to the Liouville phase, in the two other phases the central 
charge does not determine the universality class of the model. For example 
two different microscopic realizations of matter with the same large c : the 
multiple-spin model with n = 2c spin species and the Gaussian scalar model 
with D = c fields, have different susceptibility exponents 7 equal 1/3 and 1/2 
respectively. By using some general arguments one can show that there exists 
a possible series of models with positive values of 7 in the range (0, 1/2) [175 



We denote the value of the susceptibility exponent in this series by 7. The 
models in the series are related to the unitary models 7 = — 1/m, m = 2, 3 . . . 
by : 

7 = . (4.9) 

7-1 

According to the picture advocated in such surfaces with 7 look like trees 
of weakly touching bubbles which themselves describe surfaces with 7. The 
first model in this series 7 = 1/3 corresponds to random surfaces consisting of 
bubbles of pure gravity (7 = —1/2) weakly touching each other. A microscopic 



realization of such a model is the multi spin model |T9[. In this model magne 



tized domains cover the whole surface of the bubbles. Each bubble contains 



aligned spins which decouple from geometry within the bubbles [^. Effec- 
tively each bubble behaves therefore as pure gravity. The domains of aligned 
spins try to minimize the mutual contact border so that neighbouring bub- 
bles contact by a very narrow neck. Such a surface indeed looks like a tree of 
bubbles. ^ g. 

A candidate for a model with 7 = 1/4 (m = 3) being the next in the 
series should have bubbles with 7 = — 1/3 ie with the c = 1/2 matter which 
is realized by the Ising field or the Majorana fermions in continuum. The 
exponent 7 consistent with 1/4 was measured at the phase transition of the 



model of random surfaces with extrinsic curvature [^. Extrinsic curvature 



terms can be obtained by integrating out fermions from the super-string theory 



8J and therefore it is tempting to speculate that this model indeed inherits 
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some fermionic properties leading to 7 = 1/4. 

One can find a realization of the situation described above in the matrix 



model by introducing a touching term Such a touching term allows sur 



faces with a given 7 to touch each other. It has a certain coupling constant, 
x, controlling the number of such touchings. The value of the exponent 7 as 
a function of the touching coupling x, stays at the Liouville value 7 as long 
as X is smaller than a critical value Xcr- At the critical point Xcr the value of 
7 jumps to 7, and then for x > Xcr, to the branched polymer value 1/2. The 
phase with 7 is marginal in this model similarly as the marginal phase in the 
model of branched polymers. 

The phase transition at Xcr has been recently analyzed |^| in terms of 



the renormalization group flow in the parameter space There are two 

fixed points for c < 1 : one associated with the Liouville phase and the other 
with the branched polymer phase and there is a multi-critical point related to 
them on the a; = line. At c = 1, the two fixed points merge at a; = 0, and 
for c > 1 disappear from the real plane and become complex conjugate. 

When c grows they slowly depart from the real plane. Renormalization group 
trajectories in the real plane that pass in the vicinity of the complex points 
look very much like in the limiting c = 1 case, where 7 = 0. Only if one starts 
a renormalization group trajectory at /i very close to the critical value /icr 



lAi - Z^crl ~ exp 



(4.10) 



one can avoid the critical slowing down from the trajectory passing close to 
these complex points. One has to go to lattices of the size N ~ l/\fj, — fieri 
to see the branched polymer value 1/2. When c is increased, the branched 
polymer regime comes closer. This explains a long standing puzzle of values 
of 7 computed at finite volumes ||1T|, The values of 7, measured 

numerically, smoothly increase with c and reach 1/2 within the error bars only 
at c around 5. The model with the touching term has been recently simulated 
numerically ||8^ confirming the scenario of p6|. 



This renormalization group picture furnishes completely our understand- 
ing of the behaviour of the system at the border line between the Liouville 
phase and the branched polymer phase. Below we will discuss the behaviour 
of the system at the border of collapsed phase. As we show, it can be under- 
stood in terms of the balls-in-boxes model. 

With each triangulation one can associate a distribution of vertex orders 
{qi}- The opposite statement is not true, since in general it can be more 
than one triangulation associated with a given distribution Denote the 

number of such triangulations by N{{qi}). The idea is to substitute the 
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sum over triangulations by a sum over the order distributions {g^j's : 

T.7^y^T.Nm)--- (4.11) 

and to approximate by the mean-field formula : 

N{{q.})^p{qi)...p{qN), (4.12) 

where p{q) are one-vertex terms. The measure term H the partition 

function ( [4 .71 ) has the same factorized form in gj's, so it in a sense enhances the 
one-vertex terms in comparison with multi-vertex terms neglected by the ap- 
proximation. The matter fields also contribute to the one-vertex terms. In the 
large D limit this contribution is approximately q^ ^^"^ . For triangulations 
with vertices and genus h 

Y^qi = 6N -2 + 2h (4.13) 

i 

as results from the Euler relation. This constrains the sum of the vertex orders. 
One obtains the following approximation of the partition function ( ^.7| ) : 

Z^Yl ^(^i) • • • PMKQi + ...qN-{QN-2 + 2h)) (4.14) 

where P{q) = p{q)q~^^'^~°'. One recognizes the balls-in-boxes model discussed 
in the previous section. Thus one expects that for large D the transition to the 
collapsed phase, where singular vertices arise, occurs at the line a ~ —D/2. 
Generally for any D, one may force the transition to the collapsed phase 
by taking a large enough. Similarly to the collapsed phase of the branched 
polymer model, the number of triangles in the nearest neighbourhood of the 
singular vertex grows with total triangulation size. Following this, the average 
distance between points does not grow with the triangulation size. Therefore 
one can say that the Hausdorff dimension is equal infinity. 

Let us discuss now the two other phases of the model. A common feature 
of the Liouville gravity and the branched polymers phase is that the typical 
surfaces in these phases are highly branched. As a corollary of the KPZ relation 
the number of surfaces with the disc topology whose boundary is a minimal 
triangular loop is : 

N{A) ~ A^-^e'^-^ , (4.15) 

where 7 is the susceptibility exponent for spherical surfaces. A triangular loop 
on the spherical triangulation with A triangles divides the surface onto two 
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discs with the triangular loop boundary : one with area a and the other with 
A — a. The number of realizations of such a situation is 191 



MA{a) N{a)N{A-a). (4.16) 

The smaller of the discs, looks as an outgrowth on the rest of the triangulation. 
Such an outgrowth is called a minbu : minimal neck baby universe [^. When 
one considers triangulations with fixed size A, one can skip terms independent 
of a in the expression ( [4. 16] ) since they give a normalization factor. What 



remams, is : 

MA{a) ~ (^a{A - a)^ ^ ~ (^a(l - aM)^ ^ . (4.17) 

In the range 1 a <^ A/2, the distribution of minbu sizes scales as : Myi(a) ~ 
a'^~'^. In this range the tree of minbus is self similar. For 7 > the average 
minbu size : / aMA(a) ~ A"^ grows with the size of the triangulation. 

There is a close relation between the branching structure of the baby 
universes and the fractal structure of the random surface. The relation has 
been investigated by means of the real space renormalization group method 



9^, The idea is to define an elementary blocking transformation of the 
renormalization group using self similarity of the tree of baby universes. This 
idea can be practically realized by cutting off the last generation of minbus and 
calculating a change of the scale associated with the rescaling of the minbu 
tree. An interesting outcome of these studies is that the change of the average 
distance between points (r) on the decimated triangulation is related to the 
change of scale : 



{r)2 [{A), 



{r)i V A 



(4.18) 



where the subscripts 1 and 2 refer to the ensembles before and after the renor- 
malization group transformation. The blocking is performed on the fixed area 
Ai ensemble. The exponent u was found to approach u = l/dn = ^/4: for large 
lattices showing the existence of a close relation between the fractal structure 
of random surfaces and the branching structure of baby universes. 

One can determine the Hausdorff dimension directly from the scaling of 



the puncture-puncture correlation functions ||21], The correlation function 
has been found analytically for pure gravity. The resulting form shows indeed 
that it is a function of one universal parameter having the form x = r 
where dn = 4:. 
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The fractal structure of random surfaces in the presence of matter is not 

There are different theoretical predictions for 



yet fully understood ||, |9 



the value of the Hausdorff dimension. The diffusion equation in the Liouville 
theory combined with the De-Witt short distance expansion of the heat kernel 
tells us g, m 



d 



H 



' 725^+ 



(4.19) 



An alternative result obtained using the Hamiltonian formalism where one 
identifies the geodesic distance with the proper time is [DH] : 



d 
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H 



1-C+ J(25-c)(l-c) 



(4.20) 



Numerically the measurements of the matter fields with the central charge c 
in the range < c < 1 suggest that the Hausdorff dimension d^ is equal to 



four irrespectively of the matter dressing . The are two ways of measuring 
the Hausdorff dimension numerically. Either one counts the number of lattice 
points n{r) at a distance r from a random point and averages it over points 



and surfaces and then one fits the result to the formula n{r) 



dH-1 



alternative way is to measure the average distance (r) between all points on 
the lattice and then determine dn from the scaling formula (r) ~ A^^'^" ( p.l2|) 
for large lattice sizes A For pure gravity c = and for c = —2, the two 
estimates give the same values [^. In general they need not be equal. The 
definition based on the average distance (r) is related to the universal scaling 
r/N^/'^H and the mass exponent ( |2.8| ) and therefore it is closer in spirit 
to the continuum physics. For the numerical purposes the scaling argument is 
usually modified by a small finite shift r— >r + 5: x=(r + a)/N^/'^^ which 
can be neglected in the large limit [BO] . Introducing the shift for finite A^'s, 



improves the fit quality in the whole range of r 

The results of very extensive simulations can be summarized as follows. 
The numerical measurements disagree with the transfer matrix prediction 
( |4.20| ). For gravity in the presence of the matter fields < c < 1 the value 
of the Hausdorff dimension seems to approximately equal four for all c < 1. 
The values predicted by ( [4.19|) lie a bit outside the error bars of the measured 
values, but contrary to this, for c = —2 one obtains a perfect numerical agree- 
ment dn = 3.574(8) with the diffusion formula ( [4.19|) p9| . This shows that 
this part of the two dimensional theory is yet weakly understood. The fractal 
structure is being currently intensively investigated. 
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Another quantity characterizing the fractal structure of random surfaces 
is the branching dimension ds- It measures the scahng of the average number 
no{r) of disconnected pieces of the ball's boundary with radius r : 

(no(r)) ~ r'^^ (4.21) 

The branching dimension was measured numerically for c = 1 and was es- 
timated to be larger that 2.5 showing indeed a big rate of branching of the 
surfac(fl g. 

The scaling dimensions (Ib and dn can be derived from the loop distri- 
bution function E^. This distribution carries the most complete information 



about the fractal structure. The loop distribution function p(r, L) is defined 
as the average number of loops of length L on the boundary of ball with radius 
r. More precisely p{r, L)dL is the average number of loops with lengths in the 
range L to L+dL. The loop distribution was found analytically for pure gravity 
by the transfer matrix method [^. For A = oo it reads : p(r, L) = 1/ B?f{x), 
where x = L/R^, and /(x) = (x"^/^ + l/2 x"-'^/^ + 14/3 x^/'^)e~^. The moments 
of the distribution : 

(L") = ^ dLVpir, L) (4.22) 

are : 



(Li) =cirVv^ (4.23) 



(L") = c^r^"- for n > 1 

where e is the short distance cut-off {eg lattice spacing) and c's are constants. 
The zeroth moment corresponds to the number of loops of the ball's boundary 
and the first moment of the distribution corresponds to the length of the 
boundary. The two scaling dimensions c?// = 4, = 3 are related to the 
singular part of the loop distribution p (^4.23|) . 

To summarize this section. The two dimensional theory is in a very good 
shape. The phase structure is determined. One understands the behaviour of 
the system at the critical lines between phases. The Liouville phase, related 
to 2d quantum gravity, can be studied by a variety of methods. For the time 
being, the only open question is the fractal structure of surfaces in the Liouville 
phase. 

^The results of Bal were obtained from fits without the shift r ^ r + 5. 
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Figure 5: Ergodic set of local transformations of dynamical triangulations 



5 Monte Carlo simulations 

A bonus from the lattice regularization is the possibility to perform Monte 
Carlo simulations. Computer simulations provide a powerful experimental 
tool to investigate nonperturbatively statistical systems. In many cases, where 
analytic techniques break down, computer simulations are the only method to 
study a model. This for instance is the case for higher dimensional random 
geometries. 

The basic idea behind the computer simulations is to implement a Mar- 
kov chain in the space of configurations with the stationary distribution pro- 
portional to . The chain is determined by the transition probability p{\ 
2) between any two configurations 1, 2. In practice, one proposes a simple 
scheme, called the elementary step of the algorithm, of modifying a current 
configuration to obtain its successor in the chain. One can show that the 
detailed-balance condition imposed on the transition probabilities : 



e 



-S{1) 



p{l ^ 2) = e"^(2)p(2 ^ 1) (5.1^ 



and the ergodicity of the elementary steps suffice for the algorithm to generate 
configurations from the stationary distribution . 

An ergodic set of local operations to simulate two dimensional dynamical 
triangulations is shown in figure ^ The flip operation preserves the area. It 



is ergodic in the ensemble of triangulations with a fixed area [41, 100]. The 
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other two operations in the figure which remove or add a point of order three, 
change the area and allow for extending simulations to the grand-canonical 
ensemble. The three moves form a set of moves called in general the {p, q) 
moves [|101||. The first argument p corresponds to the number of triangles 



before the transformation, and q after it. Note, that the four triangles, p + q = 
4, form a tetrahedron when glued together [24]. This observation was used to 
generalize these transformations to higher dimensions, as we shall see later. 
There is also another ergodic set of transformations, called the split and joint 
operations, also used in update schemes | |102| ]. The two sets are equivalent. 

The standard algorithm was extensively tested in two dimensions. The 
distribution of triangulations generated by the algorithm is in perfect agree- 
ment with the analj^ic formula for the diagram enumeration. Numerical mea- 
surements of the critical exponents of the statistical models < c < 1 are in 



excellent agreement with the KPZ results ||103|| . The agreement extends be- 
yond the critical region as is shown by the comparison of the numerical results 
for the Ising model on the dynamical triangulations and the analytic results 
of the two matrix model |7^, [L04|| . This can be treated as a proof for the 
practical ergodicity of the Monte Carlo algorithms. This practical proof is in 
a sense stronger than the mathematical proof which only states the existence 
of a Markov chain between any two configurations which may of course not 
suffice for practical purposes. 

The {p, q) moves have a natural generalization to higher dimensional sim- 
plicial manifolds [101]. In four dimensional case there are five moves p-|-g = 6, 
p = 1, . . . , 5. Geometrically they can be viewed as a substitution of p 4- 
simplices from the triangulation by q new simplices being a complementary 
part of the boundary of a 5-simplex. One can show that the {p, q) moves are 
equivalent to the Alexander transformations [ |10CI| ] known to be ergodic in the 
set of combinatorially equivalent simplicial manifolds with a fixed topology. 
As in two dimensions, one can show the practical ergodicity of the algorithm 
|105|| . One has to keep in mind, however, as prompted in [|106|| , that the rec- 



ognizability conjecture states that for some topologies a fraction of manifolds 
accessible by a Markov chain may scale to a number less than one in the 
large volume limit. This may systematically bias the analysis of the finite size 
scahng. 

The ergodicity and the detailed balance condition leave a large freedom 
for the invention of optimal algorithms. The local update schemes are known in 
general to suffer from the slowing down, decreasing the algorithmic efficiency. 
The reason lies in the random-walk nature of local changes - namely, many 
changes are undone by successive steps of the algorithm. A general strategy 
to cure the problem is to implement algorithms focusing directly on physically 
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Figure 6: Elementary step of the baby universe surgery on dynamical trian- 
gulations. 



important modes. This strategy has been frequently used in the standard field 
theoretical Monte Carlo simulations as for instance in the multi scale or cluster 



algorithms ||107|, |108 



As discussed, typical random surfaces from the branched polymer phase 
or from the Liouville phase are populated with baby universes forming self 
similar trees |^. The idea to use the tree structure in the update scheme 



leads to the baby universe surgery algorithm ||109|| . The algorithm is in fact 
very simple as depicted in figure |^. One finds a minimal neck on the surface, 
cuts the surface along the neck, removes the corresponding minbu and pastes 
it into a randomly chosen place on the rest of the surface. Reshuffling minbus 
speeds up the algorithm dynamics by decorrelating the tree branches. 

A typical quantity characterizing the efficiency is the autocorrelation 
time. It tells us, roughly speaking, how many sweeps of the algorithm are 
needed to decorrelate the measurements of a given observable. The rise of the 
integrated autocorrelation time r for large system sizes A is controlled by the 
dynamical critical exponent z : 

T^A' (5.2) 

In the table § we present the comparison of values of the dynamical exponent 
z for the standard algorithm and a hybrid of the standard algorithm with the 
minbu surgery for various observables measured in the model with one scalar 
field. The values of the exponent z get generally reduced if one supports 
the algorithm by the minbu surgery ||109|| . This improves significantly the 



algorithm dynamics. The baby universe surgery is also applied to simulate 



higher dimensional simplicial manifolds in the elongated phase pG, 33 . 
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quantity 


^local 


^hybrid 


d 


1.06(3) 


0.76(3) 




0.81(6) 


0.14(2) 


r 


1.4(1) 


0.50(3) 



Table 2: The dynamical exponent for the standard local algorithm and the 
hybrid with the baby universe surgery, for the simulations of surfaces with 
c = 1, for the averaged internal distance between points : d, the distance 
between two given points : d^y and the gyration radius : r. One sees a large 
reduction of the exponent z when one adds the surgery to the update scheme. 



Apart from the dynamical Monte Carlo techniques described above, in 
some particular cases there are the so-called static algorithms. They sample 
directly the static distribution without using an auxiliary Markov chain. 
Thus, they are free of the dynamical slowing down. Example of such an 
algorithm is a recursive sampling technique, applicable for c = 0, — 2. In these 
cases there are known analytic formulas for the diagram enumeration which 
allow one to directly construct the diagrams weighted by the static distribution 

e"^ pTi, [nq. 

The numerical algorithms for sampling of random geometries have be- 
come a well established method allowing for studying triangulations with sizes 
ranging up to million triangles or hundreds of thousands of 4-simplices in the 
four dimensional case. 



6 Simplicial gravity 

Simplicial gravity is an attempt to formulate the quantum theory of gravity. 
It is a part of a larger programme based on the assumption that one can apply 
the same set of fundamental principles as in field theory, to quantize gravity. 
Following these lines one tries to apply the Euclidean version of the Feynman 
formalism as already described for the two dimensional gravity ( |4.1|) . There are 
many conceptual problems related to the Euclidean formulation like for exam- 
ple the lack of formal conditions which would ensure that we can reconstruct 
the Minkowskian quantum gravity. This is, however, a part of a more general 
difficulty, namely that we do not know how to formulate the Minkowskian 
quantum gravity, so we do not know what exactly should be reconstructed. 
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The role of the topology is also not clear. And again, the problem is more 
general since we can not classify four dimensional topologies. Finally, there is 
a problem with the unboundedness of the action coming from the conformal 
mode. This problem, fortunately, is automatically cured by the discretization 
scheme. Having all these difficulties in mind let us formulate in the beginning 
a more modest aim, to define consistently the Feynman integral over the Rie- 
mannian structures with a fixed topology. Now the measure problem arises. 
It is much more pronounced in four dimensions than it is in two. One way to 
bypass it is to formulate the theory perturbatively in terms of the Gaussian 
measure which is well defined. The perturbation theory obtained in this way 
is however not renormalizable and the treatment fails. One can improve the 
convergence of the perturbative diagrams at large momenta by resummation 
techniques but then one encounters problems with unitarity [|1 1 If . 



One attempt of the nonperturbative formulation of quantum gravity is 
based on the 2 + e expansion P, |112| , |113|| . The situation is somewhat similar 



to the nonlinear sigma model | ]114| | where the perturbative treatment does not 



work above two dimensions but one can find a well defined nonperturbative 
fixed point in terms of the 2+e expansion. The 2+e expansion may in principle 
be applied to gravity but in this case one has two conceptual difficulties. First 
of all, e equal two is not a small parameter, and second of all, two dimensional 
integrals appearing as coefficients in the e-expansion are not able to reproduce 
the whole content of higher dimensional gravity. 

Another way to define the Feynman integrals nonperturbatively is pro- 
vided by the lattice regularization. In this context it was proposed in P, 0. 
This is a generalization of the dynamical triangulation approach from two 



15| , |T6|, and three dimensions [^, 0, |25|1 . 

A lattice formulation allows the use of the statistical ideas and techniques 
to calculate quantum amplitudes which in the statistical language correspond 
to some averages over the statistical ensembles. The standard procedure of 
the the statistical approach is well established. Namely, given a model one 
addresses the following questions : 

1. does the model posses a well defined thermodynamic limit, 

2. what is the phase structure of the model, 

3. can one define a lattice independent continuum limit. 

The last issue is related to the existence of a continuous phase transition and 
the infinite correlation length of some physical excitations which would be 
independent of the short-range details of the lattice. 
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Let us follow these lines in the presentation. We consider an ensemble 
of simplicial manifolds with a fixed topology. Four dimensional simplicial 
manifold consists of equilateral 4-simplices glued together in such a way that 
each two neighbouring simplices share a three dimensional 3-simplex. The 
neighbourhood of each vertex is homomorphic to the 4-ball. 

As discussed in the previous sections, the sum over one dimensional 
graphs could be generated by the perturbative expansion of the scalar field 
theory while the sum over two dimensional graphs by the matrix field theory. 
The extension of this idea to the tensor model to generate four dimensional 
simplicial manifolds does not work, though it might seem to be an evident 
generalization at a first glance. The reason is that the perturbation expansion 
of tensor models generates graphs with fluctuating topology. Even worse, the 
topology of diagrams fluctuates locally so that such diagrams do not corre- 



spond to manifolds |Tl|. For the time being the only method to investigate 
the sum over the ensemble of four dimensional simplicial manifolds are the 
numerical simulations combined with the standard statistical data analysis. 

One considers the grand-canonical ensemble of simplicial manifolds with 
a given topology. The partition function reads : 

/.o) = E c^e-^^'^^^-^^'^^f^' (6.1) 



The discretized Einstein-Hilbert action (|1.6| ) reproduces in the naive contin- 
uum limit the continuum counterpart (|1.3|). For convenience we substituted 
the number of triangles N2 from the formula ( |1.6|) by the number of vertices 
A^o- This can be done since the numbers Ni of i-simplices on the manifold are 
linearly related by the Euler and Dehn-Somerville relations, which leave only 
two independent NiS. In particular, for the four dimensional sphere A^o and 
N2 are related by A'2 = 2(A'o + A'4 - 2). 

One can rewrite the grand-canonical partition function as a sum : 

Z{k^, /^o) = E ^(^4, «:o)e-'^^'^^ (6.2) 

where Z{N4, kq) are the partition functions for the canonical ensembles which 
have a fixed volume A^4. One can go one step further and express Z{N4, hq) 
in terms of the state density z(A4, A'o) : 

Z{N,, Ko) = E e'^°^''^(Ar4, No) (6.3) 

N4 

The thermodynamic limit exists if the free energy density has a well defined 
large A^4 limit. This means that the function A in the formula 

log Z{N^, Ko) = N^{f{Ko) + A{N^, Ko)} (6.4) 
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must be a finite size correction A which vanishes for large volumes A'4 — oo : 

A{N,,Ko)-^0. (6.5) 

so that in the limit A^4 00, the free energy density log Z {N4, kq) / N4 is a 
function of Kq only. 

One can show that if the function A vanishes for one finite value of 
it does so for all other. The function A was studied for different values of the 
coupling kq |]5^, p6|, ETl In particular it was found that for k,q = the 



function delta scales as A ~ where x = 0.5(2) for manifolds with spherical 



and toroidal topology [Q. This is a numerical proof for the existence of the 
thermodynamic limit. For kq = the partition function Z{N4^, = 0) is 
a sum of all manifolds without an extra weight. This sum is an object of 



intensive mathematical studies |115, 116||. The goal of these studies is to prove 



the existence of the exponential bound Z{N4^, = 0) < e'^^*. Of course, this 
is equivalent to the existence of the thermodynamic limit of the model. So far 
there is no such proof and one has to rely on the numerical results. 

One can show, as we shall see later, that the pseudo-critical value Ko,cr, at 
which the system enters the generic branched polymer phase is finite, or more 
precisely, it is bounded from above by a finite value independent of A^4. On the 
other hand, it is known that the branched polymer phase has a well defined 
thermodynamic limit. Combining these two facts, one concludes that there 
exists a finite value of hq, for which the equation (|6.5| ) is fulfilled. This suffices 
to end the proof of existence of the thermodynamic limit. In our opinion this is 
the strongest numerical evidence for the existence of the thermodynamic limit. 
What is nontrivial here is that the pseudo-critical value Ko,cr does not move to 
infinity when A^4 is increased as would happen if there was no thermodynamic 
limit. 

Numerical simulations were performed for three topologies : the sphere 
and tori x x x , x S^. In all these cases the free energy 
density /(kq) was found numerically to have the same thermodynamic limit 

Let us now outline the basic facts gathered using the computer simula- 
tions about the phase structure of the model. 

Simplicial gravity has two phases, crumpled and elongated, named to re- 
flect their geometrical properties. The elongated phase corresponds essentially 
to the branched polymer phase . In this phase typical simplicial manifolds 



are populated with baby universes which form the generation trees. The sus- 
ceptibility exponent is 7 = 1/2. The puncture-puncture correlation function is 
in this phase given by the universal formula ( ^.16| ), ( |2.17D for branched poly- 



mers. The only dependence on kq enters the formula through one universal 
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Figure 7: Minbu trees for configurations taken at random from the elongated 
phase (left) and the crumpled phase (right). Each link on the tree corresponds 
to a minimal neck on the simplicial manifold, and the end points of the link 
correspond to two parts of the simplicial manifold on both sides of the minimal 
neck. The number of links emerging from a vertex of the tree corresponds to 
the number of minimal necks found directly on the part of simplicial manifold 
associated with the vertex. The trees carry only topological information about 
the connectivity of minbus. The vertical ordering seen on the picture results 
from visualization procedure and has no relevance for the tree structure. The 
tree on the left hand side has many generations and no distinguishable points. 
The tree on the right hand size is short and has one singular vertex drawn as a 
root of the tree. There are so many links emerging directly from the singular 
vertex that the visualization procedure failed to draw them as separate lines 
and instead drew a densely covered triangle. 
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Figure 8: The minbu-minbu correlation function measured on minbu trees 
in the elongated phase of simplicial gravity (dashed line), compared with the 
universal formula G{r) ~ 2are~"^ for the branched polymers (solid line) 

(ill)- 



parameter a as in the formula (p. 16 ). In the figure |^ we show the minbu 
trees constructed on simplicial manifolds picked up randomly from two differ- 
ent phases of simplicial gravity. The tree on the left hand side comes from 
the branched polymer phase. A vertex on the tree corresponds to a minbu 
and a link to the minimal neck between neighbouring minbus. One can de- 
fine a distance between vertices on such a tree as a number of links between 
them. Then one can measure a minbu-minbu correlation function analogous 
to the puncture-puncture correlation function The results fit perfectly 

the branched polymer correlation function ( |2.16| ) as shown in figure ^ This 
means that the minbu trees are indeed branched polymers. 

The branching structure of minbus collapses when the system enters the 
crumpled phase. The collapse is related to the appearance of a singular vertex 
on the minbu tree On the tree on the right hand side in figure the 

tree has only few generations and one vertex has many branches. This is the 
singular vertex. One can check that the number of branches emerging from 
the singular vertex grows with the total number of minbus. This is shown 
in figure To distinguish this minbu from the rest, one calls it the mother 
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Figure 9: The distribution of the minbu orders in the crumpled phase of 
simphcial gravity at kO = 1.0 for three different A^4 = 4000, 8000, 16000. The 
position of the peak corresponding to the singular minbu vertex shifting pro- 
portionally to N4. 



universe. The volume of the mother universe grows extensively with the total 
volume of the simplicial manifold The reason for this rapid growth is 

related to the singular vertices residing on the mother universe. There are 
two singular vertices at a distance one. They form a singular link [^]. The 
volume of the neighbourhood of the singular vertices grows proportionally to 
the total volume. The situation is analogous to collapsed branched polymers. 
Simplicial manifolds have infinite Hausdorff dimension in this phase. 

Let us finish the survey of properties of the model by describing what 
happens at the phase transition. The standard way of investigating the be- 
haviour of a model at a phase transition is to perform the finite size analysis 
of quantities related to derivatives of the free energy. For the transition driven 
by kq one investigates the second cumulant : 

c,(.o, N,) = - = (6.6) 

which has an interpretation of the heat capacity and is a measure of thermo- 
dynamic fluctuations. The cumulant is related to the integrated correlation 
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function of the curvature so its behaviour indirectly measures the signal from 



the two point function |117|| . In particular, if the transition is second order, 
this signal may be related to the occurrence of long range correlations in the 
system. 

In the thermodynamic limit, A^4 ^ oo, the fluctuations are expected 
to approach a A^4-independent value C2(/to) except at the transition point 
where the large N/^ behaviour is approximately given by the finite size scaling 
formula : 

C2(«:o, N,) ~ Nt'^^'fiiK, - K,^,r)Nl"n (6.7) 

where a and v are the standard critical indices. They are related by the Fisher 
scaling relation a = 2 — rf/^z/. The value of the product duT-' is bounded by 
the thermodynamic inequality : 1 < dni'- In the limiting case dj^i/ = 1, the 
transition is of first order. In this case the exponent a = 1. The prefactor 
in the scaling formula (|6.7| ) grows linearly with A^4. When 1 < duv < 2 and 
< a < 1, the transition is of second order. The prefactor in (|6.71 ) grows 
as a power of A^4, between zero and one. Finally, when dnv > 2 and a < 0, 
the prefactor in (|6.7] ) does not grow with A^4 and then C2 approaches a finite 
constant at the transition. 

Finite size analysis of the second cumulant shows that the large A^4 be- 
haviour is in agreement with the first order scaling |ET|. The observed linear 



rise with the volume of the maximum of the heat capacity is related to the 
latent heat. Hence one expects a double peak structure in the N0/N4 dis- 
tribution. Indeed, such a structure has been found (see figure 0) at some 



pseudo-critical value of kq ||51|, ^ 



To summarize, there are two phases separated by a first order phase 
transition. This means that there is no continuum theory associated with this 
critical point. We shall discuss the issue of the continuum theory at the end 
of the section. 

Before doing this let us come back to the thermodynamics. The scaling 
properties of the elongated phase are exactly the same as those of the generic 
branched polymers as shows the analysis of the minbu-minbu correlations. 
The analysis of the minbu trees suggests that the correspondence to branched 
polymers can be extended to the crumpled phase as well, when taken with 
some precaution Indeed, the distribution of vertex orders on the minbu 
trees has a peak departing from the rest of the distribution when the number 
of minbus grows. Compare the figure |^ and the figure |^. 

The application of the balls-in-boxes model can be extended beyond the 
effective theory for minbu trees. Namely, the model gives also a plausible ex- 
planation of the appearance of the singular vertices on the simplicial manifold. 
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Figure 10: The distribution of Nq at kq = 1.258 for = 32000 measured 
in the computer simulations. The vertical axis corresponds to the number of 
entries for a given Nq. 

If one repeats the same line of arguments as in the section on random surfaces, 
one obtains a constrained mean field approximation for the distribution of the 
vertex orders of simplicial manifolds. The partition function for the ensemble 
with Nq vertices and N4 4-simplices is approximated as follows : 

z{N,,No) = . ..p{qNo)S{qi + ■ ■ ■ g^Vo - ^N,) (6.8) 

{qi} 

and one ends up with the balls-in-boxes model with the density p = 5N4/N0. 
If one lowers A^^o, keeping A^4 fixed, the density decreases and one triggers 
the transition to the collapsed phase with singular vertices. The mean field 
approximation assumes the independence of orders of vertices as a first approx- 
imation. The approximation does not give any particular form of p{q). Some 
numerical estimates for the weights p{q) were given in [Q. In the standard 
numerical setup one uses the (kq, iV4) ensemble. In this ensemble, the average 
Nq grows with kq. Thus one expects the low density phase for large and 
the high density for small kq. Indeed this is the case. Moreover the (kq, ^^4) 
ensemble corresponds to the (k, N) ensemble of the balls-in-boxes model ( p. 121 ) 
which has a discontinuous phase transition with a double peak distribution of 
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M/N (figure As we have seen, this is what one observes in the numerical 
data for simplicial gravity, too (figure p!0D . Introducing some next order correc- 
tions to the mean-field approximation by taking into account the geometrical 
structure of four dimensional simplicial manifolds one can explain appearance 
of the singular links as well |^ . 

The phase structure of the model discussed in the present section resem- 
bles the phase structure of random surfaces above D = c = 1. The respective 
line in the (a, D) plane (see figure 1) crosses only the branched polymer phase 
and the collapsed phase, exactly as in four dimensions. We know that in two 
dimensions, the system enters the branched polymer phase when the entropy of 
spiky conformal configurations becomes dominant. The entropy is confronted 
with the measure term whose dominance would mean that the system is in 
the collapsed phase. There is no other possibility on the line above D = 1. 
To open the physical window for the gravity phase in two dimensions, one has 
to change the matter content of the theory by decreasing the central charge 
c = D below one. This is an important lesson. 

One can similarly expect that the phase structure of four dimensional 
simplicial gravity depends on the matter dressing of the theory Out- 
side the physical window the system is realized either by the branched polymer 
phase or by the collapsed phase. The question is now how to find the window 
for the gravitational phase. To start with, one can formulate a more moderate 
goal, and ask how to prevent the system from entering the collapsed or the 



branched polymer phase. This problem has been addressed recently ^ 
Again it is useful to refer to the analogy with the two dimensional case. The 
instability of the Liouville phase is caused by the entropy of spiky configu- 
rations of the conformal field 0, ^ |62[. One can now try to repeat the 
arguments to the effective action for the conformal mode in four dimensions 
59| . The coefficient standing in front of the conformal anomaly : 



depends on the matter content of the theory through the number of scalar 
fields Ns, vector fields Ny and fermions Np, coupled to gravity |p.l8|| . The 
—28 comes from the ghost sector. The contribution from gravitons Q'^rav has 
not been calculated since it strongly depends on the ultraviolet physics for 
which the perturbative treatment fails. The coefficient plays the role of 
the effective central charge in the theory and it enters the estimates of the 
free energy of spiky conformal configurations. The reasoning is analogous as 
for the Liouville action in two dimensions [^, ^ |62[. It turns out, that for 



less than a certain critical value : < Qlrit the system is dominated by 
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the entropy of spikes which means that the system is in the branched polymer 
phase^. Hence, contrary to two dimensions, one expects that the fewer degrees 
of freedom coupled to four dimensional gravity, the higher is the probability 
for the system to be in the branched polymer phase. An inspection of the 
formula (|6.9[) shows that the strongest increase of the conformal charge 



comes from vector fields. In the paper |Q, the value of Qcrit ~ Q'irav "^^s 



estimated to be of order unity. If this is true, this means that by adding 
few generations of vector fields, one should be able to prevent the branching 
induced by the entropical instability of the conformal factor. Indeed, contrary 
to the previous investigations where some other matter fields were used 



59| , |70i [7T| , which did not affect significantly the phase structure of simplicial 
gravity, the recent simulations with a varying number of vector fields have 
shown that, when three vector fields are minimally coupled to gravity, the 
model has no branched polymer phase. Instead a new phase is created, where 
the value of the susceptibility exponent is negative, similarly as in the Liouville 
phase of the two dimensional gravity. Now we can readdress the question about 
the existence of a critical point in the extended model with the matter fields 
and about the order of the phase transition. 

To summarize, statistical physics of four dimensional complexes with 
the simplest geometrical action depending on the number of points and the 
number of four simplices is well understood. It is, however, not related to the 
continuum physics and we understand why. We hope that the problem can be 
cured by extending the phase structure of the model by adding the appropriate 
matter fields. 



Summary 

We have reviewed statistical models of random lattices used gularization 
of the problem of summing over the internal random geometry of one, two and 
four dimensional objects. The degree of difficulty in solving the problem grows 
with the dimensionality of the system, as one might have naively expected. 
On the other hand, we have shown that there are some common mechanisms 
and features, like the geometrical collapse or the existence of the branched 
polymer phase, which are almost independent of the dimensionality of the 
problem. Indeed, although the collapsed phase looks slightly different in one, 
two or four dimensions, the primary feature, namely that it is related to the 
existence of singular vertices created as a surplus anomaly, is common for all 
cases. The same universality can be found in the branched polymer phase. 

^Note, that in two dimensions the analogous inequality is in the opposite direction c > 1. 
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The most interesting part of the theory, ie the physical window, where 
the discrete models can be related to the continuum physics, is dimension 
dependent. 

In two dimensions, in the Liouville phase, related to the continuum 
physics, the universal properties of the model, like the scaling dimensions, 
are entirely determined by a single parameter being the conformal charge of 
the theory. There is a discrete series of unitary models of conformal matter 
coupled minimally to gravity with the central charge in the range between zero 
and one which have a realization as statistical models on dynamical triangu- 
lations. The Liouville phase is basically a theory of the conformal factor. 

The situation is more complicated in four dimensions where the theory is 
much more complex and requires extending the analysis beyond the conformal 
sector. But already the analysis of the conformal sector imposes some restric- 
tions on the combinations of the numbers of various fields needed to avoid the 
conformal instability. One expects that there are some other conditions one 
has to impose on the number of generations of various fields, which lead to a 
specific mixture of fields for which the theory is well defined. 
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